Absence of the zero bias peak in vortex tunneling spectra of high temperature 

superconductors 
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The c-axis tunneUng matrix of high- Tc superconductors is shown to depend strongly on the in- 
plane momentum of electrons and vanish along the four nodal lines of the d^2_y2-wa,Ye energy gap. 
This anisotropic tunneling matrix suppresses completely the contribution of the most extended 
quasiparticles in the vortex core to the c-axis tunneling current and leads to a spectrum similar to 
that of a nodeless superconductor. Our results give a natural explanation of the absence of the zero 
bias peak as well as other features observed in the vortex tunneling spectra of high- Tc cuprates. 
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Recently there has been intensified attention on the 
electronic structure in the dj.2_„2-wave vortex state of 
high-Tc superconductors |p]-[l3|| . From the solution of the 
Bogoliubov-de Gennes (BdG) equation, many authors 
showed that in a d-wave superconductor there are 
no vortex core bound states and the local density of states 
(LDOS) at the core center has a broad zero energy peak. 
However, in contrast to the theoretical prediction, the 
scanning tunneling spectrum (STS), which is generally 
assumed to be proportional to the LDOS, at the center 
of a vortex core in both YBasCusOy^^ (YBCO) [§ and 
BiaSrzCaCusOs+a (BSCCO) |J| have revealed a num- 
ber of features which are totally unexpected in a pure 
d-wave superconductor: first, there is zero peak at zero 
bias; second, the coherent peaks at the superconduct- 
ing gap edges are largely suppressed; and third, as in a 
s-wave superconductor, localized vortex core states seem 
to exist at ±5.5meV in YBCO and ±7meV in BSCCO 

i- 

Many theoretical concepts, including spin-charge sep- 
aration, SO (5) symmetry, and time reversal symme- 
try breaking pairing states have been invoked to re- 
solve this discrepancy between theories and experiments 
[f7|, p^|JT^Jl^ . A simple theoretical explanation for the van- 
ishing zero bias peak is that the pairing state inside the 
core has dx2_y2 + idxy or dx2_y2 +is wave symmetry and 
therefore as in a s-wave superconductor vortex bound 
states exist in these materials However, there is no 
direct experimental evidence for the existence of such 
time-reversal symmetry broken states. In the SO (5) the- 
ory the vortex core can be either metallic or insulating 
depending on the doping level or other physical parame- 
ters of a high-Tc material This is consistent with 
the STS data of YBCO ff~^d BSCCO |,|. Under 
the framework of the U(l) gauge theory, Franz et al 
argued that the pseudogap behavior of the STS at the 
core center is a result of the charge-spin separation and 
what the experiments found is the excitation spectrum of 
paired spinons which are not affected by external fields. 

Nearly all theoretical analyses of the STS at the vortex 
core have assumed that the tunneling matrix element is 
a momentum independent constant. Under this assump- 



tion it can be shown that the c-axis tunneling differential 
conductance measured in the STS experiments is propor- 
tional to the LDOS of low lying excitations. However, for 
high-Tc cuprates, the c-axis hopping integral of electrons 
is highly anisotropic |14-16| and depends strongly on the 
in-plane momentum. In particular, it has the same nodal 
structure as the c?2.2_j^2-wave gap function in tetragonal 
cuprates like BSCCO. The tunneling matrix is propor- 
tional to the c-axis velocity of electrons, thus it should 
also depend strongly on the in-plane momentum of elec- 
trons. This anisotropic c-axis hopping integral has strong 
impact on the c-axis transport of quasiparticles in the 
superconducting state. It leads to some peculiar temper- 
ature dependence of the c-axis penetration depth [|_5 17 



as well as microwave conductivity 18 1 at low temper- 
atures and suppresses significantly the c-axis tunneling 
conductance [Q. 

In this paper, we present a theoretical analysis of the 
STS in the vortex core of high-Tc superconductors. We 
shall show that the c-axis STS is dramatically modified 
by the anisotropy of the tunneling matrix. Our work pro- 
vides a natural explanation for the absence of the zero 
bias peak in the STS at the core center. It shows that the 
pseudogap feature of the tunneling conductance at the 
core center is mainly due to the suppression of the tunnel- 
ing current of low energy quasiparticles by the anisotropic 
c-axis hopping integral. A critical test for our theory is 
to measure the tunneling conductance in a vortex core 
which is parallel to the Cu02 planes. The tunneling ma- 
trix parallel to the Cu02 planes is not sensitive to the 
pairing symmetry and thus can be taken approximately 
as a constant. In this case, a broad zero bias peak is ex- 
pected to exist in the STS of the vortex. This anisotropic 
feature of the STS is absent in other theories. 

To study the structure of a d-wave vortex along the 
c-axis, we have performed a self-consistent calculation of 
the Bogoliubov-de Gennes (BdG) equation on a square 
lattice. Before presenting the detailed numerical results, 
let us briefly discuss why there is a broad zero-energy 
peak in the LDOS of a dx2-y2-wsLve vortex and how it is 
suppressed in the scanning tunneling experiment by the 
anisotropic c-axis hopping integral. 



1 



For conventional s- wave sup erconductors, it was estab- 
lished many decades ago 120| that discrete quasiparticle 
states with a characteristic excitation energy given by 
A"^ /2£p exist in a vortex core, where A is the bulk gap 
and Ei? is the Fermi energy. Intuitively, this can be un- 
derstood by drawing an analogy to a simple quantum 
mechanical problem of a particle in a cylindrical well of 
radius ^ ~ Ui?/7rA and height A. However, in a dx2_y2- 
wave superconductor, the radius and height of the anal- 
ogous potential well depends on the polar angle on the 
Fermi surface. Along the four node directions, A van- 
ishes and ^ diverges. In this case the quasiparticle is 
extended along the node directions and there is no truly 
localized core states. Thus unlike the s-wave case, the 
DOS of quasiparticles in the vortex of dx^-y2--wave su- 
perconductors is finite even at zero energy. From the nu- 
merical solution of the BdG equation for a dj.2_j,2-wave 
superconductor, it has been further shown that there is 
a broad peak at zero energy in the LDOS of the vor- 
tex core. This peak mainly arises from the quasiparticle 
excitations along the four node lines. 

The STS is determined by the LDOS and the tunneling 
matrix elements. In high- Tc materials, since the c-axis 
tunneling matrix vanishes along the four node lines, the 
contribution of the most extended quasiparticles, around 
the gap nodes, to the tunneling current is completely 
suppressed. This is why the zero energy peak in the 
LDOS is absent in the STS of the vortex. 

The BdG equation for a 1^3.2 _j, 2 -wave superconductor 
on a lattice is given by 

and 

HaUn{i) = -t'^Un{i + S) - a*^^„(«), 

<5 

where S denotes a nearest-neighbor vector, t is the 
hopping constant, and n is the chemical potential. 
{un{i),Vn{i)) is the wave function of Bogohubov quasi- 
particles. We have ignored the coupling to the vector 
potential A in the hopping term in the limit of an ex- 
tremely type II superconductor. The effect of the gauge 
field is manifested in the winding number of the order pa- 
rameter in the calculation. The self-consistent condition 
for the gap parameter is 

A..- - I 5] [un{^)v:{^ + S) + ujt + <5)<(*)] tanh 

71 

(3) 

where g is the pair coupling constant. 
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FIG. 1. Local density of states at the center of the vortex 
at site (0, 0) (solid line) and at site (10, 10) (dotted line) 
on a 51 X 51. The small oscillations of the curves are due 
to the finite size effects. The energy is normalized by the 
maximum energy gap Aq = 0.6t. The broadening parameter 
is 7 = 0.05*. 

The LDOS is determined by the wavefunctions of 
quasiparticles: 

p{i,E)=J2pihn,E) (4) 

n 

p{i, n, E) = \un{i)\^L{E - En) + \vn{i)\^L{E + K), (5) 

where L(x) — ^ /i: (x^ +7^) and 7 an energy broaden- 
ing parameter. E) is a sum of the LDOS of the n'th 
vortex eigenstate at site i, p[i, n, E). Figure 1 shows the 
LDOS at site (0,0) (core center) and site (10,10) on a 
51 X 51 square lattice with open boundary conditions. 
The parameters used are 17 = 1.5 and /i = 0. In agree- 
ment with other numerical results we find that the 
LDOS at the core center exhibits a broad peak at w ~ 0. 
This peak drops very quickly with the distance from the 
core center. The LDOS at site (10, 10) behaves in nearly 
the same was as in an homogeneous d^^-yi-wave super- 
conductor. There is a coherent peak at the gap edge 
(|a;| ~ Ao) at (10,10). However, at the core center, this 
coherent peak is suppressed. 

The tunneling current between a normal tip and a su- 
perconductor at a bias voltage V is given by 

n.p 

ANip,e + eV)[nFie)~nFie + eV)] (6) 

where As{n,i,e) = 2'Kp(i,n,e) and Aisi{p,e) — 2TrS{e — 
£_p) are the spectral functions of the superconductor and 
normal tip, respectively, p is the momentum of electrons 
in the normal tip and (n, kz) are the quantum numbers of 
the vortex states. The tunneling matrix element Tnk^,p 
is proportional to the velocities of electrons along the 
tunneling direction on both sides of the junction |2^ ]: 

\Tk.n-p\'^Diez)\v;vlJ, (7) 
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where D{ez) is the barrier transmission coefRcient. is 
the Fermi velocity in the normal tip, which can be ap- 
proximately taken as a constant. The velocity of the n'th 
vortex state ^ is determined by the energy dispersion 
of quasiparticles along the c-axis. To the lowest order 
approximation of t 1. (defined below), it is proportional 
to the expectation value of the c-axis current operator in 
the n'th vortex state. 

For tetragonal high- Tc compounds, the c-axis energy 
dispersion of electrons has been shown to have the form 



(fc) = —t± (cosfcj, — cosky) coskz 



(8) 



where is the c-axis hopping constant. This anisotropic 
c-axis hopping integral results from the hybridization be- 
tween antibonding O 2p and unoccupied Cu 4s orbitals 
which act as intermediate states |15| . It is also a good 
approximation for YBa2Cu307_a; |l4[| . 

Eq. can be obtained from the following interlayer 
hopping Hamiltonian 



- h.c. 



where m is the index of Cu02 planes and Dg = 1 or — 1 
if i5 = ±x or ±y. In a vortex system, since the trans- 
lational symmetry is broken and {k,j:,ky) are no longer 
good quantum numbers, it is more convenient to use this 
real space form of the c-axis hopping integral. The c-axis 
current operator is now defined by 



Jc 



E 

n.S.S' 



h.c. 



From the expectation value of Jc in the vortex state, we 
-U14|sinfc,|, (9) 



find to be 



where 



is the c-axis velocity factor of the n'th vortex state. 

Substituting Eqs. (|^) and into (^) and integrating 
out /c^, we find that the tunneling differential conduc- 
tance is 



g{i,V)^J2^nP{i,n,eV). 



(10) 



If Vn does not depend on n, such as in an ordinary 
superconductor, g(i,V) is simply proportional to the 
LDOS p{i, eV). However, as discussed below, in high- 
cuprates Vn varies strongly with n. In this case g{i,V) is 
completely different to p{i,eV). 
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FIG. 2. The c-axis velocity factor Vn as a function of 
energy. 

Figure 2 shows the energy dependence of Vn for the 
single vortex system discussed previously. The vari- 
ance of Vn is highly non-monotonic. At low energies, 
Vn is very small. It reaches a maximum when En is 
slightly above the maximum energy gap Aq. At high 
energies, Vn drops with increasing energy. The be- 
havior of Vn in different energy regions can be under- 
stood from the behavior of the expectation value of 
(cosfca; — cosfcy)^ in each vortex eigenstate. At low en- 
ergies, Vn is small since it arises mainly from excitations 
around the gap nodes. In this energy region, Vn has 
roughly the same energy dependence as (cos k^ — cos ky) 
and its local maximum increases approximately with . 
Vn becomes large when En is approximately equal to 
Aq. This is because the maximum gap corresponds 
to the contribution of the quasiparticle around (tt, 0) 
where {cos kx — cos ky) is maximum. The high en- 
ergy eigenstates are mainly from the excitations around 
(0,0) and (±7r, ±7r). Vn drops at high energies since 
(cos kx — cos kx)^ — 4 sin^ (kx + ky) /2 sin^ {kx — ky) /2 
and sin^ {kx + ky) /2 is zero at (0, 0) and (±7r, ±7r). 

Figure 3 shows the tunneling differential conductance 
for the system studied above. The conductance at 
(10, 10) behaves in qualitatively the same way as for the 
LDOS, but its coherent peak at the gap edge is sup- 
pressed. At site (0,0) (the core center), the conductance 
differs completely from the LDOS. In particular, three 
distinct features appear in the tunneling conductance. 
(1) The c-axis velocity Vn suppresses completely the zero 
energy peak of LDOS and leads to a small tunneling con- 
ductance at zero bias. The zero bias conductance is finite 
in the figure. This is due to the broadening of the spec- 
trum by the parameter 7. In an infinite lattice system 
where 7 can be asymptotically set to zero, the value of 
the zero bias conductance should become zero since the 
tunneling matrix vanishes at the gap nodes. (2) There 
are two weak peaks at low bias. The positions of these 
weak peaks depend on Aq. For the case shown in Figure 
2, they are located at V ^ ±0.25Ao. These peaks result 
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from the competition between the decreasing density of 
states and the increasing velocity factor (on average) at 
low energies. Unlike the irregular oscillations caused by 
the finite size effects, the positions of these peaks do not 
change as the lattice size is increased. This shows that 
these weak peaks are not due to the finite size effects. The 
peak value is very small compared with the conductance 
at the gap edge \V\ ^ Aq, but it slightly increases with 
increasing lattice size. The existence of these low bias 
peaks with the absence of the zero bias peak is reminis- 
cent of the tunneling spectrum of a conventional s-wave 
superconductor, although no vortex core bound states 
exist in this case. (3) There is no coherent peak at the 
gap edge. This is different than in a conventional s-wave 
superconductor. All these three features agree well with 
the experimental observations ||^-||] . 
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FIG. 3. Tunneling differential conductance at the core 
center at site (0,0) (solid line) and at site (10,10) (dotted 
line) for the same system as shown in Figure 1. 

In the above discussion, particle-hole symmetry has 
been implicitly assumed since the chemical potential fi 
has been set to zero. For a particle-hole asymmetric sys- 
tem (p 7^ 0), we find that the differential conductance 
becomes asymmetric with respect to zero bias. However, 
the other features of the tunneling spectrum seen in the 
particle-hole symmetric system are unchanged. There- 
fore, we believe the asymmetry of the STS observed in 
the experiments is mainly due to the breaking of particle- 
hole symmetry 

In conclusion, we have presented a theoretical analysis 
of the vortex tunneling spectra of high- Tc superconduc- 
tors. Our work indicates that it is important to include 
the anisotropy of the c-axis hopping integral in the anal- 
ysis of the STS in high- Tc oxides. The tunneling matrix 
elements parallel to Cu02 planes do not have the same 
anisotropy as the c-axis case. Thus STS for a vortex per- 
pendicular to the c-axis is expected to behave very dif- 
ferently. In particular, the zero energy peak of the LDOS 
at the core center should appear in the STS in this case. 



Thus by measuring the STS parallel to Cu02 planes, we 
can obtain a better understanding of the c-axis tunneling 
experiments. 
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